Distillability sudden death in qutrit-qutrit systems under global decoherence 
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Recently Song etal, Phys. Rev. A 80, 012331 (2009), have discovered that certain two-qutrit 
entangled states interacting with multi-local decoherence undergo distillability sudden death whereas 
their locally equivalent states do not exhibit such behavior. We have found that their basic equation 
of motion missed some terms and therefore the subsequent calculations and conclusions are not 
correct /valid except that distillability sudden death may occur. We provide the missing terms and 
point out the corrections which must be done in the literature to remove confusion. Moreover, 
we have studied this feature for certain entangled states under global, collective and multi-local 
decoherence. We have found that entanglement sudden death and distillability sudden death may 
happen under global and multi-local noises. Moreover, local unitary transformations can not avoid 
distillability sudden death. 
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I. INTRODUCTION 

Quantum entanglement is a vital resource for the field 
of quantum information sciences [l[ . In a noisy environ- 
ment, it is important to consider possible degradation 
due to decoherence of any initially prepared entangle- 
ment. Recently, Yu and Eberly [2, y] investigated the 
dynamics of two-qubit entanglement under various types 
of decoherence. They found that the decay of a single 
qubit coherence can be slower than the decay of qubit 
entanglement. The abrupt disappearance of entangle- 
ment at a finite time was named "entanglement sudden 
death" (ESD). Clearly, such finite-time disentanglement 
can seriously affect applications of entangled states in 
quantum information processing. The initial report of 
ESD in two-qubit entangled states in a specific model 
was later explored in wider contexts and in higher di- 
mensions of Hilbert space when the qubits are replaced 
by qutrits or qudits [|, i, i, 0, @, ©, M, E| • The exper- 
imental evidences for this effect have been reported for 
optical setups 12 1 and atomic ensembles [l3j . 

For qubit-qubit and qubit-qutrit systems Peres- 
Horodecki criterion divides the set of quantum states 
into entangled and separable states [3, [HI]. However 
for higher dimensional bipartite systems such character- 
ization is not easy [l6j]. Nevertheless, bipartite entan- 
gled states can be divided into free- entangled states and 
bound-entangled states [H, [13] • Free-entangled states 
can be distilled under local operations and classical com- 
munication (LOCC) whereas bound-entangled states can 
not be distilled to pure-state entanglement no matter how 
much copies are available. Bound entanglement itself can 
not be used for quantum information processing, how- 
ever it may activate teleportation fidelity [l^]. Bound 
entanglement manifests the irreversibility in asymptotic 
manipulation of entanglement [19( . Bound-entanglement 
constructed by purely mathematical arguments have ex- 
istence in physical processes [20j | . 

Recently, it was proposed that free-entangled states 
may be converted to bound-entangled states in the pres- 
ence of local decoherence [2~jj . In particular, it was shown 



that certain free entangled states of qutrit-qutrit systems 
become non-distillable in a finite time under the influence 
of local decoherence. Such behavior has been named as 
distillability sudden deathfDSD). We have found that 
the calculations done in [21[ are not correct. Neverthe- 
less, the main idea is still working, i. e., there are certain 
free-entangled quantum states which dynamically con- 
vert to bound-entangled states. We have rectified the 
equation of motion for the physical system and pointed 
out the possible cause of confusion. In addition, we have 
studied this problem for more general decoherence, i.e., 
when the qutrit-qutrit system is subjected to both col- 
lective and multi-local noises. We have found that dis- 
tillability sudden death may occur under multi-local and 
global decoherence. However, if qutrit-qutrit states are 
subjected only to collective noise then there are certain 
states which never undergo distillability sudden death 
and entanglement sudden death. 

This paper is organized as follows. In Sec. [Til we dis- 
cuss the physical model and the basic equation of motion 
along with its solution. We point out the corrections to 
be made in the published literature. We discuss briefly 
the idea of distillability sudden death. In Sec. IIII1 we 
demonstrate the possibility of distillability sudden death 
for both global decoherence and multi-local decoherence. 
Finally, we conclude our work in Sec. IIV1 



II. 



DYNAMICS OF QUTRIT-QUTRIT SYSTEM 
UNDER GLOBAL DECOHERENCE 



Our physical model consists of two qutrits (two-three 
level atoms for example) A and B that are coupled to a 
noisy environment both singly and collectively as shown 
in Figure [JJ The Hamiltonian of the system (h = 1) can 
be written as 0] 

H(t) = -| [b A {t) + b B (t) of + B(t)(a? + of)] ,(1) 

where [i is gyromagnetic ratio and u z denotes the de- 
phasing operator for qutrits A and B. The stochastic 
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FIG. 1: Two qutrits A and B are collectively interacting with 
the stochastic magnetic field B(t) and separately interacting 
with the stochastic magnetic fields b A (t) and 



magnetic fields refer to statistically independent classical 
Markov processes satisfying the conditions 

(6i(t)&i(f)> = ^6(t-t'),(B(t)B(t')) = ^5(t-t'), 



(h(t)) = 0, (B(t)) = 0, 



(2) 



with (• • • ) as ensemble time average, Ti and T 2 denote 
the phase-damping rates for multi-local and collective de- 
coherence, respectively. 

Let |2), |1), and |0) be the first excited state, sec- 
ond excited, and ground state of the qutrit, respectively. 
We choose the basis {|2,2>, |2,1), |2,0), |1,2), |1,1), 
1 1, 0) , |0,2), |0, 1), |0,0)}. The time-dependent den- 
sity matrix for two-qutrit system is obtained by tak- 



ing ensemble average over the three noise fields, i.e., 
p(t) = («p rt (t)}», where p st (t) = U{t)p{Q)W{t) and 
U(t) = cxp[— i J Q dt' H(t')\. The dynamics of density ma- 
trix can be given by operator sum representation [|| @] 
as p(t) = J^j Kj(t)p(0)Kj(t), where Kj are Kraus oper- 
ators [22| that preserve the positivity and unit trace con- 
ditions, i.e., KjKj = The most general solution 
of p(t) interacting with global (multi-local + collective) 
decoherence under the assumption that the system is not 
initially correlated with any of the three environments is 



given as 



P(t) 
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i,j=l k=l 



AB^ F p E A ])m{E A F B D AB )m 



where the terms describing the interaction with the 
local magnetic fields b A (t) and bs{t) are = 
diag(l, j A (t), lA {t)) ® I 3 , E$ = diag(0, u A (t), 0) ® 
Ha, E£ = diag(0, 0, u A (t)) <g> I 3 , Ff = I 3 O 
diag(l, 7B (<), 7b(*)), #f = Us ® diag(0, w B (t), 0), and 
F 3 B = I 3 O diag(0, 0, u B {t)). The ter ms describing 
interaction with collective magnetic field B(t) involve 
the operators D^ B = diag(7(t), 1, 1, 1, 7(f) , 1, 1, 1, 7_(f)), 



AB 



diag( Wl (i), 0, 0, 0, w 2 (t) , 0, 0, 0, wa(t)), D^ 



AS 



D 

diag(0,0,0,0,w 3 (£),0,0,0,w 3 (t)). The time dependent 
parameters are defined as 7a = 7b = e~ Flt / 2 , uj A (t) — 

Vl-ii(t) , U> B (t) = y/l-ij( t), 7 ^C - r ^ 2 , Wl (t) = 

v/l-7 2 W, ^2 (Q = - 7 2 Wy/l-7 2 W, and w 3 (f) = 

(1-7 2 W)V1 + 7 2 W- 

Eq. ([3]) contains 27 operators of the form, Gi = 
Si Fi Di,G 2 = Ei F x D a> . . ., G 27 = £ 3 F 3 D 3 , and in- 
deed they satisfy the condition J2 n =i G n = \- The 
matrix form of Eq. ([3]) is given as 
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We note that if multi- local fields are turned off, i. e., ja = 
^b = 1, and qutrit-qutrit system is subjected to collective 
noise, then decoherence free subspaces (DFS) 0, 0] do 
appear in this system. Similarly, we can concentrate on 
the dynamics only for multi-local decoherence by turning 
off the collective magnetic field, i. c., 7 = 1 in Eq. 

We point out here that Eq. (1) of 21] is not correct 
as the authors studied the case of multi-local noise 



(with 7 = 1) and defined the general solution as 

P(t) = Eij=i F ]( t ) E l( t )p(°) E i F 3 (compare it with 
Eq. ©). This summation contain four terms instead 
of nine. They missed the matrices £3 and F3. In fact 
the equation (p(t) = £ 2 J=1 Fj{t)E}(t)p(Q)E i F j ) was 
written in [6(, where the authors studied the existence 
of entanglement sudden death for qudit-qudit (d <g> d) 
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systems. In that case, the solution must be p(t) = 
J2tj=i F j( t ) E l( t )p(0)E l F J and the corresponding 
operators must be E\{t) — diag(l, 7a, 7a, ■ ■ ■ ,7a) ® Id, 
E 2 {t) = diag(0,WA,0,...,0) <g> I d ,...,E d (t) = 
diag(0,0, 0, . . . ,oj a ) ® l d , and F 1 (t) 
I d (g> diag(l,7s,7s, ■ • ■ ,7b), ^(t) = Id ® 
diag(0, us, 0, . . . , 0),. . = I d ® diag(0, 0, 0, . . . , u B ). 

It is quite straight forward to show that with Ei and 
Fj defined in [a], the partial transpose of isotropic 
states always contain one negative eigenvalue and the 
conclusions that isotropic states always exhibit entangle- 
ment sudden death would not be correct. However, the 
conclusions of Ref. Q are correct and we assume that 
Ann-Jaeger only wrote the wrong matrices but worked 
with correct ones. The proof of this assumption can be 
found in Ref. Q, where again Ann- Jaeger wrote wrong 
operators in their Section 2, however they provide the 
correct matrix as Eq. (2) of Ref. 0], which validates their 
conclusions. However, such confusion should be rectified 
for broader interest of readers. In fact, the matrices Ei, 
E2, and E3 for a qutrit were used by Ann- Jaeger in 
their study of entanglement sudden death in qubit-qutrit 
systems @. Song etal, [2l[ probably have taken the 
results of Ref. [|| and discovered the very interesting 
feature of distillability sudden death for qutrit-qutrit 
systems. As a result of this confusion, Eq. (3) of Ref. [2l[ 
is wrong. However, as we demonstrate in Sec. IIIII that 
some of their conclusions still survive, in particular, 
free-entangled states undergo distillability sudden death. 
However, the local unitary transformation can not avoid 
distillability sudden death. 



III. DISTILLABILITY SUDDEN DEATH 

Before demonstrating the possibility of distillabil- 
ity sudden death, we discuss the characterization of 
bound-entangled states. It was proved [l7| that bound- 
entangled states have positive partial transpose (PPT) 
and hence non-distillable under LOCC. There is not 
a unique criterion to detect bound-entangled states. 
Even for qutrit-qutrit system, there are various bound- 
entangled states and a single criterion is not capable 
to detect all of them (sec for example [23| and refer- 
ences therein). However, one can use the realignment 
criterion [24[ to detect certain bound-entangled states. 
The realignment of a given density matrix is obtained as 
(p R )ij,ki = Pikji- For a separable state p, realignment 
criterion implies that \\p R \\ < 1. For a PPT-state, the 



1 can prove the 



positive value of the quantity \\p 
bound-entangled state. 

Let us consider a particular initial state given as 
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(5) 



where 2 < 

1/3(|0,1)(0,1| 

1/3(|1,0)(1,0| 



a 



< 5. the 
|1,2)(1,2| + 
|2,1)(2,1| + 



separable states 04 
|2,0)(2,0|), and <r_ 
|0,2)(0,2|), and 



l/v^GO, 0) + I LI) + |2, 2)) a maximally entangled state. 
It was shown [18j that p a is separable for 2 < a < 3, 
bound-entangled for 3 < a < 4, and free-entangled for 
4 < a < 5. It was claimed "pi")] that the partial trans- 
pose of time-dependent Eq. ([5]) in the range 4 < a < 5 
always has a negative eigenvalue in finite time. We have 
already mentioned that this claim is a result of an in- 
correct equation of motion. The fact is that the possible 
negative eigenvalues of partial transpose of Eq. ([5]) are 
given as 

Ai 

A 2 

A3 



167 8 7 A7 |+4a 2 -20a + 25 , 



167 8 7a7b +4a 2 - 20a + 25 



1 

42 



167 A7 | + 4a 2 -20a + 25) . (6) 



These eigenvalues obviously become positive at a finite 
time. It can be observed that if we turn off the multi- 
local fields and turn on only collective interaction, then 
these states always has a negative value, which implies no 
entanglement sudden death and no distillability sudden 
death. However, in the presence of global decoherence 
or multi-local decoherence, entanglement sudden death 
and distillability sudden death may happen as we demon- 
strate below. 

The negativity [25| for this state can be easily calcu- 
lated as 

N(p a (t)) = max[0, -Ai] + max[0, -A 2 ] + max[0, -A3] .(7) 

The negativity N(p a (t)) is plotted for various values of 
a in Figure [2j where we have taken Ti = T 2 = T for 
ease. Figure[5]shows that except for a = 5, the negativity 
becomes zero at finite times leading to a PPT state which 
might be entangled. The state with a = 5 becomes PPT 
only at infinity, which is definitely a separable state. 
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FIG. 2: The negativity is plotted against the decay parameter 
Ft for a — 4.3, a = 4.6, and a = 5. 

Let us now demonstrate that the time evolved den- 
sity matrix p a (t) do undergo distillability sudden death. 



This could be most easily demonstrated by checking en- 
tanglement of p a (t) after it becomes PPT. This relation 
is plotted in Figure which shows that under global de- 
coherence an initial free-entangled state becomes bound- 
entangled at a time Ft^ ~ 0.1422. The entangle- 
ment of PPT-state is verified by the positive value of 
\\p R {t)\\ - 1 in the range 0.1422 < Ft < 0.1764. How- 
ever, the realignment criterion fails to detect the pos- 
sible entanglement after time Ttu ~ 0.1764. To con- 
clude the separability or entanglement after this time, 
one could extract three " 2 (g) 2" density matrices spanned 
by { |1, 1), |1,0), |0, 1), |0,0> }, { |2, 2), |2, 1), |0, 2), |0, 1) }, 
and { |2, 2) J2.0), |1, 2), |1, 0) } from p a (t) and check their 
properties [211 ]. 
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FIG. 4: The negativity and the realignment criterion is plot- 
ted against the decay parameter Ft for a — 4.3. The state be- 
comes bound-entangled at time FtN ~ 0.1422 and it remains 
bound-entangled at least in the range 0.1422 < Ft < 0.3437. 
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FIG. 3: The negativity N(p) and the realignment criterion 
\\p R || — 1 is plotted against the decay parameter Ft for a = 4.3. 
The negativity is zero at time Ttjv ~ 0.1422, nevertheless, the 
state is bound-entangled in the range 0.1422 < Ft < 0.1764. 

We could turn off the collective field, i.e., 7 = 1 
and observe the distillability sudden death under multi- 
local decoherence. Figure 2] shows that p a (t) loses its 
negativity at a finite time rtjv ~ 0.1422 and becomes 
bound-entangled as the realignment criterion has a pos- 
itive value for the interval 0.1422 < Ft < 0.3437. Hence 
p a it) experiences distillability sudden death under multi- 
local decoherence as well. It is surprising to note that 
the dynamics of negativity is not much affected by turn- 
ing off the collective field but the realignment criterion 
is considerably affected. We could imagine that intro- 
ducing additional noise cause faster decay of entangle- 
ment, nevertheless, in this case non-distillable entangle- 
ment showed much respect to additional noise compared 
with distillable entanglement. 

Let us now consider the locally equivalent state <r a , 
which may be obtained from Eq. (0 with local unitary 
operation U = I3 0, with 6 = |0)(1| + |1)(0| + |2)(2|. 
Then we have 



U Pa U^ 



a 



7 



where the separable states are 
|1,2)(1,2| + |2,1)(2,1|), and a_ 
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1/3(|0,0)(0,0| 
1/3(11,1X1,11 



(8) 



|2,0)(2,0| + |0,2)(0,2|),and |#+) = 1/V3(|0, 1) + 11, 0) + 
|2,2)) a maximally entangled state. The local unitary 
transformations do not effect the static entanglement, 
however they may have a profound influence on the fu- 
ture trajectory of entanglement [lfj. The possible neg- 
ative eigenvalues of the time evolved partially transpose 
matrix are given as 



m = ^(^5- v /4a 2 -20a + 25 + 16 7 l7|J-(9) 

Again we observe that in the presence of only collective 
interaction, 773 is always negative which implies no entan- 
glement sudden death and no distillability sudden death. 
Whereas in the presence of either global decoherence or 
multi-local decoherence, the entanglement sudden death 
and/or distillability sudden death may happen. 

In Figure O we have plotted the negativity and the re- 
alignment criterion against the decay parameter Ft. We 
observe that distillability sudden death occurs in this case 
as well. The states becomes PPT at a time Ft^ w 0.0948. 
Nevertheless, the states are bound-entangled in the in- 
terval 0.0948 <Ft< 0.2686. The state a a {t) is fragile 
compared with p a {t) in the sense that a a (t) becomes 
bound-entangled at earlier times, i.e., distillability sud- 
den death comes earlier as compared with dynamics of 
p a (t) (compare with Figure [3]). However, in this case the 
realignment criterion detects bound-entangled states for 
a broader range. 

FigureOshows the negativity N(a a (t)) and ||cr^ (t) |j — 1 
plotted against the decay parameter Ft for a specific value 
of parameter a. We observe that a a (t) becomes bound- 
entangled at Ft ps 0.1422. Figure[6]is exactly like Figure 
51 which means that the local unitary transformation has 
no effect on the dynamics of negativity and dynamics of 
bound-entangled states for this particular case. 
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FIG. 5: The negativity N(a) and the realignment criterion 
||<r fl || — 1 is plotted against the decay parameter Ft for a = 4.3. 
The state becomes PPT at an earlier time TtN ~ 0.0948, nev- 
ertheless, it is entangled in the range 0.0948 < Tt < 0.2686. 
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FIG. 6: The negativity N(cr) and the realignment criterion 
||cr [| — 1 is plotted against the decay parameter Ft for a = 4.3. 
The state becomes bound-entangled at time TtN ~ 0.1422, 
where negativity has zero value. 

Finally, let us consider the isotropic states defined as 

/0p = p|* + )(* + |+(l^][ 9l (10) 

where < p < 1 . The isotropic states have the property 
that their PPT region is always separable (see Ref. [23[ 



and references therein). Therefore, we realize that these 
states do not suffer distiilability sudden death, however 
they may suffer entanglement sudden death. The possi- 
ble negative eigenvalues of time evolved partially trans- 
pose matrix are given as 

£1 = -P- 3p 7 4 7A7s) , 

6 = g (l -P~ ip 1 A 1a1b) , 

& = ~(l-p-3p7i7l) • (11) 
We can see that isotropic states for < p < 1, do not 
suffer entanglement sudden death under collective deco- 
herence. However in the presence of global decoherence 
or at least for multi-local noise, they always exhibit en- 
tanglement sudden death for < p < 1. 

IV. SUMMARY 

We have revisited the entanglement dynamics of bi- 
partite systems under global, collective and multi-local 
noises. In particular, we have studied the distiilability 
sudden death of qutrit-qutrit systems. We have found 
that certain free-entangled qutrit-qutrit states becomes 
bound-entangled as a consequence of purely dynamical 
process. We have found that distiilability sudden death 
may occur under the combined action of multi-local and 
collective noises. Distiilability sudden death can also re- 
sult under multi-local decoherence. We have observed 
that local unitary transformations can not avoid this 
phenomenon. Therefore, we conclude that qutrit-qutrit 
states evolve into three main types of dynamics: (i) They 
loose their entanglement only at infinity and hence never 
suffer entanglement sudden death and distiilability sud- 
den death, (ii) They loose their entanglement at a finite 
time but never undergo distiilability sudden death, (iii) 
They first undergo distiilability sudden death and then 
sudden death of entanglement. 
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